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Abstract 

The aim of this paper is to give a result concerning the stability properties of the solutions 
of magnetohydrodynamics equations at small but finite Reynolds numbers. These solutions 
are found using the or-effect: this method gives us solutions which are highly oscillating 
spatially on the scale of the underlying flow but are growing on a larger scale depending on 
a parameter e. We prove instability results for a dense subset of initial velocity field of the 
flow at given Reynolds number. 
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1 Introduction 

Dynamo theory first interested physicists almost one century ago, when Sir Joseph Larmor [fT4ll 
wondered how a magnetic field could be generated by a rotating body like the Sun or the Earth. 
Since then, many results have come, both in mathematics and physics, to understand this kind of 
phenomena. 

The physical model behind them is a coupled system derived from the interaction of a conduc- 
tive fluid, which follows Navier-Stokes equations, and an electromagnetic field, which follows 
Maxwell equations and Ohm's law. The interaction of the two physic quantities comes from the 
fact that a magnetic field has an effect on a moving conductive fluid, which creates in return a 
magnetic field. This kind of system is known as magnetohydrodynamics (MHD). 

There are mainly two outcomes for MHD systems: the magnetic energy can either collapse 
(exponentially), or stay/grow. In the latter case we say we have a dynamo, and it is the effect 
we observe in bodies like the Sun. However, the understanding of dynamo mechanisms is not 
simple. Indeed, the most "intuitive" way we imagine the conductive fluid in rotation and the 
magnetic field (that is a 2D rotating fluid in the kernel, and a magnetic field that looks like the 
one of a static magnet) is actually a well known anti-dynamo situation. Anti-dynamo theorems 
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have been widely studied, and some summary of the knowledge in the domain has been made in 

im. 

The kind of situation we want to study there is that of the second kind, that is dynamo situ- 
ation (and more precisely growth). Positive results began to emerge only later in the 50's, even 
though it was already well known that it is what occured in the Sun. In the 90's M.M. Vishik [|23l 
proved that fast dynamos (i.e. exponential growth) require Lagrangian chaos in the fluid. Those 
results spread the idea that the velocity field has to be chaotic enough to expect a dynamo effect. 
Because of this "negative" result, different methods were tried to find instability results, linked 
to singularities: Y.B. Ponomarenko [17J studied concentration phenomena to get dynamo result 
in the case where Laplace force is neglected (i.e. imposed velocity field) and his results were 
extended by Gilbert fTOl and successfully experimentaly checked [|6l. More realistic case of the 
Ponomarenko dynamo was studied by D. Gerard- Varet and F. Rousset AH those results give 
only asymptotic results. 

In this paper we will focus on another singular dynamo mechanism, called alpha effect, first 
introduced by E.N. Parker [|T6l , which is based on scale separation: we formally decompose the 
magnetic and the velocity fields into two parts: a fluctuating part, evolving on a turbulent scale, 
and a mean part, evolving on a larger scale. The mean value is broadly the spatial average over a 
sphere of radius much larger than the typical distance on which the field is turbulent, but much 
smaller thant the distance on which we want to observe the dynamo effect (in what follows, these 
distances will be 1 and l/s, so scale separation will take place at small s). Basically, this gives a 
decomposition of the velocity and magnetic fields under the form 

u = u{x, x/s, t) + u{x, t), 
b = b{x, x/s, t) + b{x, t), 

where u and b have null mean value on the second variable. The idea behind alpha effect is that 
the averaged induction term V a (m a 5) (we will specify those notations right below) that will 
appear in the MHD equations will permit us to obtain a destabilizing effect on the magnetic field. 

The alpha effect was first formally studied by G.O. Roberts [fT9l [TSll . in the framework of 
periodic flows. His results were discussed by H.K. Moffatt ifTSl . S. Childress [|2l and A.M. 
Soward [|20l that considered a larger class of velocity fields and range of parameters. This effect 
has been since experimentaly confirmed by R. Stieglitz and U. Miiller [|2TI . 

The aim of the present paper is to provide a rigorous justification of this mechanism, in a 
nonlinear framework. We will improve substantially former results of D. Gerard- Varet [171 [H. 

Before stating more precisely our results, we will first recall the system we will work with: 

The generation of magnetic field B by an electrically conducting fluid of velocity U and 
magnetic diffusivity v„, follows the equations 

{ d,B^V a{U aB)- v^AB = 0, 

I V • 5 = 0. ^ ^ 

These equations are derived from the Maxwell's equation and Ohm's law (see second section of 
lEl, or [[TTI for more detailed derivation). Since other relevant physics quantities that appear in 
the Maxwell's equations can be obtained from B, it is enough to work with only that. 
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Furthermore, the incompressible Navier-Stokes equations that drive the evolution of the con- 
ducting fluid under the force created by the magnetic field write 

p^,u + pU■vu = -Vp + r^^u + j-^(yAB)AB + f, 

V • M = 0, ^ ^ 

where p is the pressure field, / is an additional forcing term, p is the fluid mass density and 77 
is the fluid shear viscosity (see [3 J again, keeping in mind the incompressible condition). The 
additional forcing term will be chosen to generate some steady velocity field (and thus will be 
independent of time), see below. 

We study there oscillating field U of amplitude V, oscillating on a space-scale L, and with 
magnetic difi'usivity v^. After proper nondimensionalization of time, space, pressure and mag- 
netic field, we then have the full incompressible MHD system: 

r dtB-V Aiu aB)- ^AB = 0, 

' ^ ' R„, ' 

dtU +U -VU + Vp~ j-AU = (V A 5) A 5 + /, 

V ■ 5 = 0, ' ^ ^ 

V • f/ = 0, 

where Rm = ^ is the magnetic Reynolds number, and Re = ^ is the hydrodynamic Reynolds 
number. 

Here we want to find a growing solution on a larger scale L/s (where s > 0), using the 



alpha effect introduced earlier. We will search this solution as a perturbation y^] of a solution 
' of this system (and the forcing term / is chosen accordingly so that Us is a solution of the 





Navier-Stokes equation): 

U = u + Us, 



and the system in terms of b, u and U^ becomes 

dtu + u- VUs + Us ■ Vm - j-Au + Vp = {V A b) A b ~ u - Vu, 
dtb-V A iUs Ab) - V A (u Ab) - ^Ab = 0, 

V b = 0, ^ ^ 

V ■u = 0. 

The space variables for this system is denoted 9, and corresponds to the fluctuation scale L, 
1 after non-dimensionalization. In particular, we assume that 

[/, = C/,(0) eC"(T^) (7) 

and has zero-mean. 

Later on, we shall introduce a larger scale variable x = s6 corresponding to the large scale 
L/s of the growing field. 
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Our main result is an instability theorem in the Liapunov sense for system To state our 



result, we introduce the set 9i^(T^) with his usual distance: 



max (1, II t/ — ^y'll^O 
d{U,U')=J] ^^^^ (8) 

k=0 

which makes it a Frechet-space. 

For T = {Tu T2, Tg) e (R* )\ we denote 

"K^ = "H'iR/TiZ X R/T2Z X R/T^Z) (9) 
the Sobolev space of order s on the torus 

R/TiZ X R/T2Z X R/T3Z. (10) 
We have the following theorem: 
Theorem 1 Let s > 5/2 be a real number. There exists a dense subset 0^ of the unit ball of 
'K^(T^) and > such that for all Us e ^ and Rm < 7?° , for all Re > 0, the solution ^ 



of the (MHD) system is nonlinearly unstable in the following sense: 
There exists T = {T\,T2,T^) e (N*)^ and initial values 







(11) 



and Co > such that for all 6 > Q, the solution of ^ with initial value ^ ] satisfies 



Z ) (^) 



> Co (12) 



T 



for some time t^. Furthermore, we have the estimate on tg as 6 ^ 0: 

In 5 

h -6^0 -— (13) 

As expected, this theorem expresses the instability of the MHD system, up to consider a larger 
periodic box [0, Ti] x [0, T2] x [0, T3]. Let us point out that (Tj, Ta, T3) e (N*)^ so that t/, is 
also periodic on this larger scale. 
In particular, for any 

;«) . (14) 

and any 6 > 0, system ^ has a unique solution 

sC%[0,T%H^,) (15) 
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for some T* > with initial data 6 (uq, bo) . Moreover, for 5 > small enough, it is well-known 
that one can take T* = +oo. 

Up to our knowledge, our theorem is the first full justification of the a-effect, in the nonlinear 
context. 

In particular, it extends substentially the article by David Gerard- Varet. Recast in our 
variables, this article considers the case of large magnetic diffusion, that is i?^ = e « 1- Briefly, 
it is shown that for all m e N, there are solutions of © (on a larger box depending on s) that go 
from amplitude e™ initially to an amplitude rj = J]{m) > independent of s (but depending on 
m). 

Note that, as t] depends on m, this result does not yield Liapunov instability at fixed small s. 
Moreover, the study in [|7l corresponds to vanishing magnetic Reynolds number. In our theorem, 
the threshold is of order 1, which is much more satisfactory from a physical point of view. 

The outline of the paper is as follow: 

The first and main part of the paper is the constructino of an exact unstable eigenmode for 
the linearization of It involves sharp spectral analysis arguments. Then, on the basis of this 
linear analysis, the proof of nonlinear instability is given in section 3. 

2 Linear instability 

In this section we will establish linear instability for the induction part of the MHD system. In 
other words, we will exhibit an unstable eigenmode of the linear operator 

b^V a{U Ab)-^Ab (16) 

for a large class of U. This will be possible, through an accurate mathematical analysis of the 
alpha effect described in the introduction. 
Let us write down the induction equation: 

dtb-V a(U Ab)- —Ab = 0, 

^ ^ R,n (17) 

V ■b = 0, 

where U s 'H^(T^) is independent of time. 

In this section we prove the following linear instability result: 

Theorem 2 Let s ^ be a real number. There exists a dense subset 3^ of the unit ball of 
'K^(T^) and a positive real R^^ such that for all positive R,,, < R^^, and all U s 0^, there exists 
^ e (2;rQ*)^ and a number > such that for all < s < sq, there exists a solution B of (flTl) 
of the form 

b{t,e) = e''^-'e''''b'{e) (18) 
where b^ e 'H*(T^), with exponential growth when ? — > +oo (that is !R(/l^) > Oj. 
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2.1 Linearized equation 

Let U e 'K|^(T^) be a spatially periodic velocity field in (of period 1). 

Since we will look for the growth of the magnetic field on a larger space-scale 1/e, we will 
directly look for a solution b under the form 

b{e, t) = B%se, e, st) = B%se, st) + b^sO, e, st), (i9) 

where B'^ = B^{x, 6, t) depends both on the original variable 6 and a slow variable x = s6. One 
can decompose B^ as the sum of a large-scale field B^ which doesn't depend on the smaller 
variable, and a small-scale zero-mean-value field B^, i.e. B^{x, 6,t)A6 = 0, as suggested in 
the introduction. 

In view of (|4l), the new field B^ should satisfy (dropping the prime over time variable t): 

d,B' - V, A ([/ A 5^) - iV, A (t/ A 5^) - f-^,B'^ ^AeB' - I: Sti ^.v,^,^' = 0, 
V, - B^ + ^Vg B^ = 0. 

(20) 

The idea here is to find an exact solution of this new equation in {t, x, 6) (or more precisely to 
prove that there exists one that fits our needs). We will find this solution as a Fourier mode in x: 

B^ = e'^-'B^'{t,e). (21) 

Dropping the prime again, we get: 

d,B' - a{U A B') iV, A ([/ A B') + f-fB^ ^AgB' - j^i^^U ^i^e.B' = 0' 
■ 5^ + B' = 0. 

(22) 

We denote by A'^ = A" + eA^ -I- s^A^ the operator associated to the first equation of (|22l) . that is 

A« = V, A ([/ A •) + 

2 

A' = A {U A ■) + —iV^^de, (23) 

A' = -^eu, 

where A'' : •K'+2(T^) "K'^T^), such that the first equation in ^ writes 

dtB"" = -A^5^ (24) 

£ 

We want to find a solution in power series of s of that equation of type 

B%e,t) = e^'''Bl{e). (25) 

Thus, the equation writes: 

A'Bl = sA^B^O' (26) 
which is an eigenvalue problem. We will split this equation in powers of s to obtain 3 equations. 
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2.2 Study of the first term in tlie development 

The first equation, corresponding to e = 0, writes 

A^Bl = 0, (27) 

that is j-AqBq + Ve a {U a 5°) = 0, which depends only on the independent variable 6. 
We will look for solution of type 8^(6) = + 5^(6) as proposed before: 

^AeB° + V, A ([/ a5°) = -Ve A ([/ a5°). (28) 

Since the Laplacian 

Ae : 7^+2 (T^ ) c (T^ ) ^ (T^ ) (29) 
is an unbounded invertible operator, one can introduce 

A-^Ve A (f/ A •) : -T^UT^') ^ 'H^+^T^) ^ •^o(T^) (30) 

which is a compact operator of domain 'Kq (T^). 

Thus, by classical result on the spectrum of compact operators, the operator 

Lfi,„ = -^A, • +V, a{U A-) = Ae{^ld+A,'Vo a (t/ a •)) (31) 

is invertible as an operator from 'K()^^(T^) to 'Kg (T^) but for some countable set Cl of which 
can be written as a sequence of numbers growing to +oo. In particular it is invertible for R,„ 
smaller than a given value R^^. 

Since -Ve a {U a fijj) g 9{^{T^) for any 5° e R\ the equation (|28]) admits exactly one 






solution for every 5!] e R^. Then for any R,„ out of Cl, the equation 



L«„ 5° = -V, a{U A 5|]) (32) 

is equivalent to 

5° = ^{e)Bl. (33) 
where ^{6) is the operator from to 'Ko+^(T^) defined by 

^{6) = -L^ye A {U A ■) (34) 

That means that the operator has a kernel of dimension 3, A^: "H'+^iT^) c 'H'{T^) 
'H'{T^), that is ker A" = {B + ^{e)B, B e R^}. 

Using Lax-Milgram theorem and elliptic regularity, we know that the operator 

— A + Ve A ([/A •)-/iId=A°-iId (35) 

is invertible for large enough A if seen as an operator from "Hq"^^ to TYq. Thus the resolvant of 
A° is compact in at least one A, and all spectral projections are of finite dimension, that is A° 
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has eigenvalues of finite dimension, isolated and they are the only spectral values of (see [fT3]| 
p.187). 

Let cr be a smooth path surrounding and no other eigenvalue of A". For s sufficiently small, 
since A'^ varies continuously on s, cr doesn't cross the spectrum of A"^, and the operator 

P' ■= ^ ( {A' - Ald)-^ dA. (36) 

is the spectral projector associated to the eigenvalues of A"^ which tend to as e tends to 0. (0 
is an eigenvalue of degree 3 in the case of A", but this eigenvalue will vary - and maybe split - 
around for positive e, without crossing cr if e is small enough). Note that the path depends on 
the choice of ^. 

Morover, still for e small enough, /"^ also varies analytically on s. 

Furthermore, by definition, is the spectral projection of A" on its eigenvalue 0, and com- 
mutes with A" (lini again, thm 6.17 p 178). 

On ImP", A° has then only as an eigenvalue. Since ImP^ is of finite dimension, A'^ is 
nilpotent on this same space. Until the end of the subsection, we define = A*'| inipo. 

Let n be the integer such that 

kerjTl" c ker{Jiy c • • • c ker{jf)" = ImP", (37) 

We will now prove that n = \, that is ker = ker(j?l°)^: 
Let B E ker(j?l°)^. Then we have the equivalences: 

{Jl'^fB = « ^^i^B) + 'Wb ) = 0, 

« = B' + ^{e)B' and Wb = 0, (38) 
« ^B = and Wb = 0, 
« Jl'^B = 0. 

Thus ker = ker{Jl^f = ImP". 

2.3 Simplification of the equation 

Back to (|26l) . we look for solutions of that equation in ImP^, that is in the sum of the proper- 
spaces corresponding to the spectral values of A'^ surrounded by cr: 

Developping P^ as the beginning of a series in e, P^ = P^ + sP^ + s^P'^ , we have 

pej^sps _ pO^OpO ^ gplji^OpO ^ gpO^Opl ^ gpO^lpO ^ ^2^/^^ (4Q) 

where both P'"^ and A"^ depend analytically on s. 
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Since and A° commutes and A"/'" = (see results of previous subsection), the first three 
terms vanish, thus leaving us with the equation: 

In the next subsections, we look for a solution of this equation for e = 0. We shall then 
conclude using the continuity of eigenvalues as functions of s. 

2.4 Case e = 

The first order in e of the equation (|4TI) is equivalent to the system 



A^B = AB 
B = ^{e)B 



(42) 



where = a {U a ■) — Xlf^i ^t^er The mean part of the first equation writes 



a{U aB) = AB, (43) 

that is 

a{U A ^{e)B) + a{U aB)= AB, (44) 

^ y ' 

=0 

A (aB) = AB, (45) 
where a is a constant matrix (that depends on U and R,„), 

ab= \ U{e) A ^{e)b, (46) 

which we will study in the next subsection. 

Thus we now have to find the eigenvalues of the operator a {a-), that is to find the eigen- 
values of the matrix A^a, where 

-1^3 

A^ = I -/^i I . (47) 



-i^2 i^i 



2.5 Study of the matrix a 

The a matrix is defined for any e R.^ by 



ab= \ U{e) A ^{e)b = [ U{e) a fcVe A {U{e) A b)) . (48) 

We shall write this matrix as the sum of his antisymmetric part or'* and symmetric part 
(which are both real). 
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Now we compute the eigenvalues of those matrices: 
Since is antisymmetric, his eigenvalues are 0, +i[|y;7[| where Vu is a vector of M? such that 

ybGR\a^b= Vu Ab. 

For , we will apply a slight variation to U so that the matrix changes to one of distinct 
non-zero eigenvalues (this doesn't change what we said about the antisymmetric part). 



2.5.1 Symmetry of a 

We have in the expression of a: 



- L-ye A{UAb)=(^-^ Id ( - A-^ V, a{U A b)) = (J- Id -^^ j^b, (49) 

where 

^ = -A-'Ve A (t/ A •) : K^' ^ K^'- (50) 

For small enough (such that 1 /R,„ is bigger than the highest eigenvalue of the compact oper- 
ator £/, say R,„ < 27?° where > - we will justify this factor 2 as a precaution at the end of 
the subsection), we can express the resolvent as a series: 



+00 



- L-ye A ([/ A Z7) = 2 R':;'^"+'b = Rl^"b. (51) 

n=0 «=1 

Furthermore, (l50l) gives in Fourier series, foxk i^Q and / e 'H'^ (T^): 

= -^\k\-\^ik) A Uik') A f\k k')^ . (52) 



Note that, by definition of \ a7(0) = 0. 
thus 



U{6) A y{R^£^fb 

+00 « 



(53) 



J/A£/"i(0) 

Recursively, we can compute U a £/"b(0) and 

«=1 i;;i,'A:;=0 \ I "I \ \ I 11 / / / 



(54) 
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(56) 



where the second sum is actually restricted to tuples where all partial sums m, = Xly^i ^; 

non-zero (otherwise there is a term of the form £/''b{0) which is zero by definition of £/). 
This last equality writes: 

00 

ab = Y,R>^"^'^b, (55) 

n=l 

where 

T!'=ikj=o \ I « II \ \ I II 

and a^"^ have the same symmetry as n (by recurrence, or see [[T]|). 
2.5.2 Eigenvalues of a 

The first symmetric matrix in power series is a^^^ : 

Let U E 7-{'f^{T^), and 77 > 0. We want U such that a^^\U) is of distinct non-zero eigenvalues 
and d{U, U) < 77, where d is the usual distance induced by semi-norms on ^^(T^): we define 
now 

Uj{e) = 2 U{k)e^"''\ (57) 

For large enough j, we have d{U^, U) < 77/2. Denote then U = + where (5,- > is 

determined later, and 

'sin(2;r(j + 1)63) + cos(2;r(j + l)e2)\ / sm{2n{j + 2)^3) 

V = I cos(2;r(j + 1)^3) , = i sm{2n{j + 2)0i) + cos(2;r(7 + 2)^3) 
sin(2;r(j+ 1)02) / \ cos(2;r(j + 2)^1) 
cos{2n{j + 3)62) 

V = I sin(2;r(j + 3)ei) 1. (58) 



sin(2;r(j + 3)^2) + cos(2;r(j + 3)6* 



Then: 



.J f -men 1/ J Oif A:i ?t 0orA:2 ^ 

V {k)y = J e Vy{e)de = I ^i,-2,v.3.3cos(2;r(7+ l)03)d03 otherwise ^^^^ 

_ J if (^3 J + 1 and ^3 it - j - I) or ki or k2 
\ \ otherwise 



Thus 
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and similarily for the other ones: 

jSki=0^k2=oi^k3 = -j-2 — Sk3=j+2) 
^ (^) ^ I 2^ki=oSk2=o{^k3=j+2 + Sk3=-j-2) + f^fe=0<5fe=o(<^yti=-y-2 " ^^1=7+2) 

2^k2=oSk3=o{^ki=j+2 + ^ki = -j-2) 



^Sk2=oSk3=o{^ki = -j-3 — Ski=j+3) 



(61) 



(62) 



For y (^) to be non-zero, we need that kj = and ki e {0, j + /, —j — i} when / i. 

Then ^^(-A:) a (i^ a (V^(A:) a 5)) = as soon as A;i ?t or A:2 0. If both are null, then 

V\k) = but if A:3 e {0,7 + 1, - 1}, and V^{k) = but if A:3 e {0,7 + 2, -7 - 2}. Thus the 
product is always zero. 

Only remains now non-crossed terms: 

^\-k) A {i^^ A {y\k) A B)^ = V\-k) A V\k). {i^^ ■ . (63) 
For it to be non-zero, we need that ki = 0. If furthermore k2 i= and ^3 0, then (k) = 0. 



If its =0andA:2 ^0 




Thus for A:2 = 7+1, V\-k) a ^^(A:) = i 



If its ^0andA:2 = 

|((^*:3 = -7-l - 5^:3=7+1)' 
V'{k)= I i(43=,-+l+43 = -,--l) 



Thus for A:3 =7+1, V\-k) a V\k) = \\ 



and thus 




(64) 



12 



We do the same for the other ones: 



k*0 



7 + 2 



), Y,n-k)-nk)(i^,-B 



Finally, for all i e 11,3]], since the spectrum of U-' and V are disjoint: 



and then 



0,(2) (c7) =a^^)^uJ)- 



(2), 



Si 



1 



'B. 
(65) 

(66) 



^0 0\ ^ /I 0\ ^ /I 0^ 
1 --^ -^ 1 

j + ^ \o I j + ^ \o I j + ^ \0 0. 



(67) 



where the 5, can be chosen as small as we want. Thus a^^\U) is as, near as we want from a matrix 
of distinct non-zero eigenvalues, and it's the same for (U) for R,„ < R'^^ (NB: the set in which 
Rm can be taken will vary continuously with U, thus the precaution we took at the beginning of 
the subsection to be sufficiently far from the limit) 



2.5.3 End of the proof 

In the sequel of this subsection, we now work with U instead of U, droping the tilda. Thus, 
has distinct non-zero eigenvalues. In Fourier transform, we now have to find the eigenvalues of 
the matrix A^a^ + A^a^ , where 

-ft ft 

= / I ft ^1 I . (68) 



-ft ft 



Let there he P e such that 




'Pa^P = a2 \ =D, (69) 



and write ^ = 'P^. The matrix A^D has Aq = eigenvalue, and two opposite eigenvalues given 
by 



A+ = + Y ^fa2as + ^20^3^1 + ^i^iai- (70) 

Those eigenvalues are real for any ^ if all eigenvalues of have the same sign. Otherwise, there 
exists a cone with vertex outside of which they are. Furthermore, for ^ in the complement Q. 
of this cone, X+ stays positive, and X_ stays negative. Thus, the eigenvectors of A'^D depend 
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regularly of ^. From now on, we fix some C, in this good set £l. We can also assume that 
^ = P^G(2;rQ*)3. 

Let be an eigenvector associated to Then A^Dfi = a {D0) = A+/3, so ^ ■ /3 = 0. Thus, 
we have = ^ a (7 a /3) = • y)/3, where 7 is a vector such that VZ? e R\ 'Pa^Pb = 

y A b. Thus /3 is also an eigenvector of A^'Pa^P, associated with the eigenvalue — ■ 7). 

Now changing back to ^ = P^, we have that 

= 'PA^P, (71) 

so that /3 is an eigenvector of 'PA^a^P + 'PA^PD associated to the eigenvalue A+ — ■ 7) = 
/l+ — • Py). And thus P/3 is an eigenvector of A^a^ + A^a^ associated to an eigenvalue of 
positive real part, that is a growing mode solution of the first term of (|4TI) . 

2.6 Conclusion to the linear case 

We have proven now the existence of an eigenvector b*^ of the equation 

for £ = with positive real part, that is an exponentially-growing mode of the form 

B\x, e, t) = bl{e, ^)e(-^+-'-f r)m-f-x (73^ 

of the induction equation for given velocity field Us and small enough R,„ (and e = 0). Since 
the operator ^p^A^P^ varies analytically on s for small e, and has distinct eigenvalues for e = 0, 
so do eigenvalues and eigenvectors of this operator, and thus the result remains valid when e 
is positive sufliciently small (at fixed ^), and we have now an eigenvector bo of the induction 
equation for small enough s. 

This concludes the proof of theorem |2l 

Remark 1 Most of what we said stays valid if we take Rm out of a discrete set of values instead 
of only "small enough". Only the fact that a{U) is as near as we want to a matrix of distinct 
non-zero eigenvalues remains a problem in that case. 

Pondering on theorem |2l we will prove in the next section the instability of the full nonlinear 
(MHD) system ©. 

3 Nonlinear instability 

We will now prove nonlinear instability of the system dH). For that we will adapt a method 
developped by S. Friedlander, W. Strauss and M. Vishik in [5] (itself adapted from Y. Guo and W. 
Strauss in [fT2]| ) to prove nonlinear instability given linear instability. The condition s > 3/2 + I 
in the main theorem is used only in this section to obtain an inequality on the nonlinear terms 
similar to that in the article. Actually, using the parabolic feature of the MHD system, it would be 



14 



possible to lower our regularity requirements (see for instance Friedlander et al BUl in the context 
of the Navier-Stokes equations). We stick here to regular data for simplicity of exposure. 

Let Us e ^ belonging to the dense subset of linearly unstable flows described in the previous 
section. 

We introduce Leray operator P to get rid of p. We rewrite system Q as 



where 



+ L,^ = 2(^,'^), (74) 
^ = ' ^^^^ 

is the linearized operator around {^q^^ ^^'^ 

is the nonlinearity. 

We shall consider this system in Hj, T to be specified later. We also denote the space 
on the torus. Since we will work in these spaces until the end of the proof, we will drop the T in 
the notation. We have the following a priori estimates for the nonlinear term: 

Every term of every component of Q{'^ ,'^) can be written as a product of one component 
of ^ and one of V'^ (composed by P). Thus, we have the first estimate 

||2(^,^)[|l2 < C||^ • V^[|i2 < C[|^||L2||V^[|ix,. (78) 

For r = \{s + n/l + 1) > n/2 + 1 (here n = 3), we have that 

IIV^IIloo < W'^Ww, (79) 

and 

W'^Ww <C\\^\\l,\\'W\\l^J' (80) 

for T] such that r = (1 — ri)s (i.e. rj = \ — ^ = \ — j;)- 
Thus, we have: 

112(^,^)1^2 < cil^ll^+^ll^liy (81) 

We shall prove the following theorem, which is a mere reformulation of the first theorem: 

Theorems Suppose s > 3/2 + 1, small enough and Rg > 0. Then there exists a vector 
T = {T\,T2,T-i) with integer coordinates such that the system (I74l) is nonlinearly unstable in 'Hj 
in the following sense: 

There exists a growing mode = ^'^^ H'^oH'K^ = 1 of the linearized system and a 

constant Co > depending only on s,n,p such that for all 6 > 0, the perturbation of the 
linear solution d'^o satisfies \ ^{ts)\'n^ > Cq for some time tg. Furthermore, tg -^s^q 
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Remark 2 This theorem only proves that the solution reaches high values infinite time however 
small the initial condition is. In particular it doesn 't say wether it explodes in finite time or is 
defined for all time, since both situations are possible in the result of the theorem. Actually, for 
small enough initial data ( that is small enough 6 ), it is classical that the solution of the system 
doesn 't explode in finite time and thus the approximation of tg for vanishing 5 is valid. More 
details in / |221/ (Theorem 3.8) 

Assume in contrary that the system is nonlinearly stable while spectrally unstable, that is that 
for any growing mode of initial value = ( ) with \\'^q\\'U' = 1 in any box T of the system 

(i.e. an eigenvector of associated to an eigenvalue with positive real part) and for any £ > 0, 
there exists a 5 > such that for all t > 0, the solution (t) of the system with initial values 
6% satisfies < £• 

3.1 Proof of the nonlinear instability 

Let bL{6) = e^^''^b^{6) be the unstable eigenmode given by Theorem 2. We fix e < eo in such a 
way that 



belongs to N^. This is possible as ^ e (2:?rQ* )^ (take s to be the inverse of a large integer). In 




(82) 




p = max {%A, A e Spectr(Lv)} . 



(83) 




Moreover, by standard properties of the spectral radius of e^\ we know that 



lim lie' 



.'^^ II 1/' 



(84) 



Thus, for all ?7 > 0, there exists a Q such that for all f > 0, 



(85) 



Then, the solution of (|74l ) with initial data 5Uq can be written 



(86) 
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where ^ satisfies 



= 

Using Duhamel's formula, we can write the solution under the form: 



Jo 



Define 



Ts = sup I? > 0, < t, \\'^{s)\\l2 < i(JeP-^ll^ollL2 j , 

where A is the coeflacient in the growing linear mode. Then, since we supposed that ^ 
remains bounded, and ||^/|| ~ 6ef"\\'^o\\, we have that Tg < oo. 
And for all t < Tg, 

Using a priori estimates stated earlier, we have for 77 = ^ — ^ = 5 ~ ^- 

Jo 

Jo 

< C„5i+''e'(i+^ V-"— (e'5P _ 1) 
W 

Thus, getting back to ^ : 
By definition of Tg, we have: 

that is 



2 

2' 
^5 >ln 



Denote by tg the quantity 



?5 = In 



2C'5''ei-''(l +77)7 77p' 
17 



Thus ts < Ts, and at this time ts 



Thus II'^IIlz (and H^H-kO can be bounded below by a value independent of 6, which contradicts 
the asumption, and proves the theorem. 

3.2 Estimation of the growing time 

We will now estimate the time in which ||^ ||<h» becomes of order 1. Denote by Tq this time: 

To = sup{t > 0, V* < t, \\'^{s)\\w < 1}. (99) 
We use again the definition of T^: 

Ts = supjr > 0,V5 < t, \\'^{s)y < ]^6eP'\\'^oy]. (100) 
This time, Ts might be infinite. For big enough A as in the former proof, we have again that 

T, > ,„ 1 (101, 

Denote again by tg the right-hand side: 

' \2C'6''{l+T])J T]p ^ ' 

Thus in time = infjr^, Tq} and with the constant 

i„f{l. ^11*.,,. (J^) '■''}, (103) 

we have 

ll^(?o)b>Co, (104) 

and 

\2C'{l + T])J Tjp p p 
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